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1. Introduction 

Entanglement is a central aspect of modern discussions of quantum mechanics. 
Recognized in the early days[^, it was exclusively a matter of foundations for more 
than half a century [ 2 . Its status changed with the birth of quantum information 
theory[|3J, when it become a resource to be understood and utilized in tasks like 
teleportation[|S] as well as in some algorithms of quantum computation [|HJ [7]. 

One natural task is to look for entanglement sources. The most usual strategies 
are to manipulate individual qubits and to control their interactions in order to 
generate entanglement, or to post select a special part of a quantum state, using 
our knowledge of some property to create entanglement in other degrees of free- 
dom. Other tasks include the search for the consequences of entanglement for the 
properties of matter 

In this article we will revise a proposed strategy to obtain spin entangled fermion 
pairs from a degenerated Fermi gas. We will address the natural difficulties in 
applying such ideas to electronic Fermi gases and propose the use of stored neutrons 
in order to experimentally investigate entanglement in Fermi gases. 

2. Pair entanglement in Fermi gases 

In this section we review two works, refs. IIP) , which are the central contribution 
we presented in CEWQO. Details can be obtained in the cited references. 
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In Ref. [03 the conditional state of a localized pair of fermions from a degener- 
ated Fermi gas was studied. Consider |$o) as the ground state of the Fermi gas. In 
a second quantization notation, one has 

i*o>=n b ^)i°>> a) 

S , p 

where h\ (p) is the creation operator of one fermion of spin s and momentum p, and 
the product is over all modes with momentum smaller than the Fermi momentum p f . 
If one considers the two body density operator (up to normalization), conditioned 
to the detection of a pair of fermions at positions r and r' , then[ITT] 

pf^tv = < $ ol *l (r') *1 (r) *V (r') * s (r) |$ > , (2) 

where 

r Pf d 3 v 

*.(»•)=/ 7-^e^b.(p) (3) 
Jo (2,it) 

is the annihilation operator localized at r. 

The operator can be considered as a two qubit density operator (see Ap- 
pendix of Ref. [EH) an d application of Peres-Horodecki criterion [ 112) shows the 
existence of an entanglement distance: if two fermions are extracted from the Fermi 
gas separated by less than this entanglement distance, then their spins are entan- 
gled. One physical point must be stressed: the position detection must be spin 
non-destructive for the scheme to make sense. One limiting case is very intuitive. 
If both fermions were detected at the same position, as they only have position and 
spin degrees of freedom, their spin state is necessarilly the singlet state \][) — |J,|), 
which is maximally entangled. In Ref. [E3 it is shown how the increasing in the 
fermions separation corresponds to more mixing with a classically correlated state, 
which makes entanglement to vanish at and above this entanglement distance. En- 
tanglement distance is essentially given by kj , where kf is the Fermi wavenumber. 

One step forward was given in Ref. [ 110] when the perfect localization condi- 
tion was dropped out, and the notion of fuzzy mesurements was introduced. Fuzzy 
measurements are characterized by a distribution D (r — r") and detection (annihi- 
lation) operators 

*- w = /|^/|F' r " D(r - r '' )b - w - (4) 

Using this notion, entanglement was quantified and its dependence on the separa- 
tion and on the broadening of gaussian detectors was studied (note that eq. © 
corresponds to the situation when D is a delta function). However, no pratical 
situation was devised. 

It must be stressed that fuzzy measurements must be seen as coherent detectors, 
in the sense that they have specific annihilation operators associated to them. They 
are not a simple bunch of incoherent localized detectors. Microscopic detectors tend 
to be closer to fuzzy measurements, while macroscopic detectors tend to add up 
incoherent signals. 
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3. Electronic Fermi sea 

The simplest fermions one can think about are electrons and the simplest situation 
for them to be confined in a non-interacting Fermi gas is the case of conduction 
electrons in a metal. 

From a theoretic and abstract viewpoint, it can be done. The difficulty in 
proposing a real fuzzy measurement for such a gas comes from the entanglement 
distance they have, of the order of a few angstroms [Ej- One should then be able to 
address a pair of conducting electrons as close as such a distance, without disturbing 
its spin state. This is not a simple direct task, and some other strategy needs to be 
devised in order to extract entanglement from the Pauli principle. 

4. A proposal: Neutrons 

For a confined Fermi gas, fc^ is proportional to the density of fermions [I14|. So, the 
root of the difficulties with conducting electrons is their high density. Neutrons can 
be used to create a Fermi gas suitable for the extraction of spin entangled pairs. 
We will describe such a proposal in this section. 

The experiment we will propose keeps some similarities with an interesting 
recently performed measurement of evidences for neutron antibunching[ ITS] , In 
that experiment a beam of thermal neutrons is monochromatized by a cristal. The 
beam then incides in a beamsplitter and the coincidences on both outputs of the 
beamsplitter are recorded as a function of the optical path difference. A small deep 
is registered, which can be understood as a result of the fermionic antibunching 
(for each pair, the triplet spin component antibunches, while the singlet component 
bunches). The smallness of this deep can be interpreted as an indication that the 
vast majority of the recorded coincidences arc accidental, in the sense that the two 
recorded neutrons were not really paired in the beam (in other words, they were 
independent), even when the optical paths agree. This is natural for a reasonably 
high flux of the order of 3000 neutrons per second. 

The one and very first difference in the experiment we propose is the previous 
creation of a confined Fermi gas. This is made to keep the situation the closest pos- 
sible to the idealized Fermi gas, where boundary conditions imply the discreteness 
of the single particle levels. Ultra-cold neutrons [I16j are a good example of confined 
Fermi gas, but other alternatives are also possible. 

Now we want a physical way of realizing the detection operator (J3J. Note that 
the D distribution imply a wavepacket-like detector. This can be approximated by a 
sequence of a slit and collimators, where the emergent neutrons will have reasonably 
well defined position and momentum. The details of the slit and collimators will 
define the best model for D. A naive model is just a gaussian profile. 

Naturally, neutrons will be emitted by this "source" in a statistical way. The 
flux can be estimated as Nrce, where N is the number of trapped neutrons (decres- 
cent with time), r is the ratio between the area of the hole and the total area of the 
storing vessel, c is the rate of colisions of neutrons with the walls, and e < I is an 
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efficiency coefficient reiated to ioses in the collimation process or in the detection 
(also effects of antibunching can be included in this factor). If one stores N = I0 5 
neutrons in a bottle of one litre [ 1 1 7| . c can be avaliatted as c « 50 s _1 . For a hole 
diameter of less than 1 mm, one gets r » 10~ 5 . As e depends on experimental 
details, one can say that roughly the flux will be of the order of few neutrons per 
second. The experiment intends to work with essentially simultaneously emitted 
neutrons. Such a concept depends on the definition of a simultaneity time window. 
This has to do with both, fundamental and practical aspects. From the first point, 
the broadening in momentum distribution of the emitted neutron defines a time 
scale which can be associated to its presence. From the other one, the response 
time and accuracy of the detectors play an important role. Suppose we allow a 
coincidence time of the order of 1CP 4 s. This value is large compared to the usual 
neutron detectors. From the values above, one can consider to have a mean value 
of the order of 10 -3 neutrons per "pulse" , which implies a coincidence rate of the 
order of 10 -6 per counting or 10~ 2 coincidences per second. These numbers sugest 
the experiment to be feasible, despite difficult. The vessel should be fulled and de- 
tections made through a time like one hour (something around 30 pairs of neutrons 
would be available). Then the preparation must be repeated and data recolected. 

Up to this moment, we were not specific of which experiment should be done. 
If one just want to prove entanglement, something like a Bell experiment is enough. 
For completeness, we suggest to make a tomography of the spin state of neutron 
pairs [ 18 . The procedure is detailed in the appendix. In any case, a beamsplitter 
should be used to separate the neutrons and spin polarizing detectors must be used 
in the equal optical path condition for recording coincidences. 

4.1. Other entanglement effects on neutrons 

We must emphasize that it is not the first time that one talks about entanglement on 
neutrons. There is a debate about some experiments on neutron scattering in which 
"anomalous" experimental cross-section may [El or may not[|251 be caused by short 
time entanglement among nucleons. Also violations of Bell inequalities were verified 
with neutrons [12 1| which must be considered as entanglement on neutrons. However, 
while in the first case, the possibly entangled nucleons are subject to a strong 
environment which disentangles the system, and in the second experiment, the 
entanglement is between two degrees of freedom of the same neutron: momentum 
and spin; the experiment which we propose will be able to generate free flying spin 
entangled neutron pairs. 

5. Conclusions 

In this article we review how non-desctructive measurements can be used to ex- 
tract entangled pairs of fermions from a Fermi gas. We discuss the difficulties in 
applying such ideas to conducting electrons in normal metals. We then suggest an 
experiment with neutrons capable of exhibit such a behaviour. It can be seen as 
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a practical advance: a possible new source of entangled particles; and as a foun- 
dational interesting step, testing novel properties of fundamental constituents of 
matter. 

Appendix: Quantum state tomography 

We want to describe a simple two-qubit tomographic procedure for neutron spins. 
We will first introduce a formalism in terms of expectation values of observables. 
Then we will convert it to neutron-countings since these are the available data. 

The central idea can be understood for just one qubit. Denote by <jq the 2x2 
identity matrix and by o~j, j — 1, 2, 3, the usual Pauli matrices (by convention, latin 
indices run from 1 to 3, while greek indices include 0). Note that {ct m } is a basis 
for the (real) vector space of (2 x 2) Hermitian matrices. Moreover, <jj are traceless 
matrices. So, any density operator can be written as 



in a unique way. The three real numbers bj form the so called Block vector , and 
the positivity condition for p implies that this vetor lies inside a unit sphere with 
respect to the usual Euclidean norm (the so called Block sphere). Eq. (0) can be 
interpreted as a tomographic decomposition when one realizes that 



which follows directly from the algebraic properties of Pauli matrices. This means 
that one only needs to measure three mean values to characterize a one qubit state. 

In terms of particle counting, it means that for each j one needs to count two 
rates: n J+ and rij^ and the tomographic coefficient will be estimated by the ratio 

" . (7) 

n j+ + rij- 

Whenever possible, the best way is to use polarizing beamsplitters (PBS) and two 
independent detectors, one in each output of the PBS. Measurements for different 
j values involve the rotation of the PBS and the inclusion of quarter wave plates. 
If the PBS is not balanced, Eq. Q must be properly adapted. 

For generalizing this to two qubits, consider the matrices S M „ = a^®a v . Again, 
the set {S M „} is a basis for the (real) vector space of (4 x 4) Hermitian matrices. In 
a more compact notation, one can write any density operator as 





bj = Tr {o-jp} = (o-j) , 



(6) 




(8) 



with the coefficients given by 



a. 



Tr{S M „p} = (S M „) , 



(9) 
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which must now be interpreted. First of all, TrjSoop} = 1 and must not be ob- 
tained. The remaining 15 coefficients can be divided into three sets: a^o refers only 
to the first neutron, aoj to the second, while to correlations. Ideally, the first 
two sets could be measured without paying atention to coincidences. However, in 
the experiment we have in mind, we want to characterize the state of neutron pairs. 
So, it is necessary to use the second detector as a trigger for the relevant one. 

As in the one qubit case, the mean values must be translated into detection 
rates, now with only coincidences being registered. The coefficients can now be 
estimated as 

(10) 



n i+ n j+ 


- n l+ 


n 3~ 




n j+ 


+ rii- 


n 3- 


n l+ n 3+ 
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n l+ n 3+ - n l+ rij- + ni-n j+ - ni-n 3 - 

a 0j — , (J-^J 

ni + rij + + ni + nj- + rii-nj + + rii-rij- 

where each coefficient and aoj will be estimated three times. Ideally these 
estimations should coincide. In practice they will coincide up to experimental errors 
and a statistical treatment can be used to refine them. 

One important difficulty we want to address is that the use of PBS for two neu- 
trons can give rise to difficulties in the registering of coincidences. In this case one 
will need to make independent measurements for each two complementar situations. 

It is important to emphasize that, in a situation like the here proposed, where 
each outcome takes a relativelly long time, it can be interesting to consider better 
tomographic strategies, like the one proposed in ref. [ I22j . 
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